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GO Abstract 

d We present computations of cuprates which accurately describe the physics arising at scales of up to several lattice constants. 
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After 20 years of research the low energy physics of 
cuprates remains a mystery. It is generally accepted that a 
three band Hamiltonian with Cu's and O's orbitals should 
contain the essential physics so we can say that there is 
consensus at the length scale of the Cu-0 distance but 
there is disagreement when it comes to describe higher 
length scales or to set up a low energy model. 

Much on the modeling in the field consists on guessing the 
low energy physics. We adopt a more conservative approach 
and try to move from the safe land of the atomic scale to the 
intermediate scale of a few atomic constants. An exact solu- 
tion is not possible but a number of techniques allow to ad- 
dress the intermediate length scale physics with accuracy. 
Rather than guessing we can accurately compute what the 
model predicts at intermediate energies and length scales 
and compare in detail with experiments. This should con- 
siderably constrain the low energy physics. If our modeling 
is successful we should be able to describe spectroscopies 
from high energies up to some infrared cutoff (below which 
our ignorance sets in) and we should be able to describe 
short and intermediate length scale physics. This includes 
the ground state energy which, in a short range model, is 
essentially a short range property. 

In the last years we have carried on this program push- 
ing in some cases our "discerning cutoff" down to a few 
meV. We have based our computations on the Gutzwiller 
approximation and a time dependent extension pQ. These 
techniques have revealed to be very accurate when com- 
pared with exact diagonalization [I|2l3| and exact results 
in onefl] and infinite]!] dimensions. An accurate technique 



supplemented with an accurate model should produce re- 
sults in accord with experiment and it does. 

Within our approximations the lowest energy solutions 
consist of O centered stripes which have approximately a 
number of holes per lattice site along the stripe v ~ 0.5. 
These textures explain the behavior of the incommensu- 
rability, the chemical potential, some anomalous transport 
experiments^ and the optical conductivity as a function 
of doping^. In order to get accurate parameters we have 
relied on LDA. New calculations in much larger systems 
which allow better resolution show small deviation from 
experiment. For example the filling of the stripes with the 
parameters of Ref. [5] is v = 0.56 holes per site rather than 
v = 0.5 as found in experiment [7] ■ This can be easily under- 
stand from the fact that the stripe filling is approximately 
given by v ~ y/ J/2C with C the inverse compressibility of 
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Fig. 1. as a function of doping for vertical O centered stripes. 
We also show the result for Cu centered (CC) stripes for d = 4. The 
inset in reports the incommensurability as obtained from the present 
calculation (line) compared with experimental data from Ref. 0. 
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the stripe and J the superexchange interaction^ 5 . LDA 
parameters are well known to overestimate the value of J 
so it is not surprising that the filling of the stripe is overes- 
timated. To improve the agreement with experiment here 
we present new results with an empirical parameter set 
based on Ref. [9]. We use A = 3.3eV for Cu-0 onsite en- 
ergy splitting, t pd = 1.3eV (t pp = 0.7eV) for Cu-0 (O-O) 
hopping, Ud = 8.8 eV (U p = 6.0 eV) for onsite repulsion on 
Cu (O), U pd = 1.0 eV for Cu-0 repulsion and K pd = -0.25 
for Cu-0 direct exchange. The CuO interactions U p d and 
K p d where added by us. K p d has been recognized as essen- 
tial to yield the correct value of J[TD] and so it plays an 
unsuspected role on fixing the filling of the stripe. We have 
found that the optical conductivity and the magnetic ex- 
citations computed within the time dependent Gutzwiller 
approximation [Tj for the undoped system are in excellent 
agreement with the experiment confirming the accuracy of 
the parameters. 

Fig. [1] shows the energy per hole as a function of dop- 
ing. Results are qualitatively similar as in Ref. [5] but now 
the stripes have v ~ 1/2 even for large systems. As dis- 
cussed before^ for doping x < 1/8 the system responds to 
a change of doping by changing the periodicity d whereas 
for x > 1/8 the periodicity gets fixed at d = 4 in a large 
doping range. The inset shows the resulting magnetic in- 
commensurability compared with experiment. 

Fig. [2] shows the charge and spin density at doping 
x = 1/8. A similar charge distribution was predicted in 
Ref. [5] and found to be in excellent agreement with a 
charge sensitive probe by Abbamonte an collaborators [XT] . 
They could not determine if the stripes where Cu centered 
or O centered. More recently Davis and collaborators [121] 
have imaged glassy stripes which indeed are centered on 
O as predicted [5]. Taken together these experiments give 
us amplitude and phase information of stripes in excellent 
agreement with Ref. [5] thus showing that it is possible 
to obtain and even predict realistic information on the 
intermediate scale physics of cuprates. 

The Fourier transform of the charge and spin distribu- 
tion determines Bragg peak weights (oc TOq) in scattering 
experiments (for the definitions see Ref. [13]) as shown in 
Fig. Disregarding the difference in cross section for dif- 
ferent processes (magnetic neutron scattering vs. X-ray or 
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Fig. 2. Open circles (p orbitals) represent Cu (O) sites. The numbers 
and the size of symbols represent the excess charge whereas the spin 
density is proportional to the length of arrows. 
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Fig. 3. Bragg weights m,Q , Uq at reciprocal lattice vectors (indicated 
by the dots) for x = 1/8 and d = 4 stripes. The polygons around each 
point are the reduced magnetic Brillouin zones. We use reciprocal 
lattice units of the fundamental Brillouin zone. Ox (Oy) are the 
oxygens in the bonds perpendicular (parallel) to the stripe. 

nuclear neutron scattering) we see that Bragg weights are 
practically 3 orders of magnitude smaller for charge than 
for spin. This is due to very soft charge distribution shown 
in Fig. [2] with respect to the spin distribution and explains 
why it has been so hard to detect charge ordering. 

Optical excitation on top of these textures and magnetic 
excitations on top of similar textures in the one band Hub- 
bard model have been found to be in excellent agreement 
with experiment confirming the validity of the intermedi- 
ate scale physics found [T^TTS] . 

To conclude our computations are based on a mean-field 
like approach and therefore can not address subtle issues as 
the absence or presence of long range order, specially in the 
delicate magnetic channel, but we can reliably determine 
the behavior at intermediate length scales improving our 
understanding of these fascinating materials. Our results 
provide strong constrains on low energy theories. 



References 

[I] G. Seibold, J. Lorenzana, Phys. Rev. Lett. 86 (2001) 2605. 

[2] G. Seibold, F. Becca, J. Lorenzana, Phys. Rev. B 67 (2003) 
085108. 

[3] G. Seibold, F. Becca, P. Rubin, J. Lorenzana, Phys. Rev. B 69 
(2004) 155113. 

[4] F. Giinther, G. Seibold, J. Lorenzana, Phys. Rev. Lett, in press 

( |cond-mat /061 1583) . 
[5] J. Lorenzana, G. Seibold, Phys. Rev. Lett. 89 (2002) 136401. 
[6] J. Lorenzana, G. Seibold, Phys. Rev. Lett. 90 (2003) 066404. 
[7] K. Yamada,ei al Phys. Rev. B 57 (1998) 6165. 
[8] G. Seibold, J. Lorenzana, Phys. Rev. B 69 (2004) 134513. 
[9] H. Eskcs, L. H. Tjeng, G. A. Sawatzky, Phys. Rev. B 42 (1990) 

288. 

[10] E. B. Stechel, D. R. Jennison, Phys. Rev. B 38 (1988) 4632. 

[II] P. Abbamonte et al, Nature Phys. 1 (2005) 155. 
[12] Y. Kohsaka et al Science 315 (5817) (2007) 1380. 

[13] J. Lorenzana, G. Seibold, R. Coldea, Phys. Rev. B 72 (2005) 
224511. 

[14] G. Seibold, J. Lorenzana, Phys. Rev. Lett. 94 (2005) 107006. 



[15] G. Seibold, J. Lorenzana, Phys. Rev. B 73 (14) (2006) 144515. 



